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Abstract 

For Z bb, we calculate all the two-loop top dependent Feynman graphs, 
which have mixed QCD and electroweak contributions that are not factorizable. 
For evaluating the graphs, without resorting to a mass expansion, we apply 
a two-loop extension of the one-loop Passarino-Veltman reduction. This is 
an analytic-numerical method, which first converts all diagrams into a set of 
ten standard scalar functions, and then integrates them numerically over the 
remaining Feynman parameters, with rapid convergence and high accuracy. We 
discuss the treatment of infrared singularities within our methods. We do not 
resort to unitarity cuts of two-point functions for calculating decay rates; these 
are useful only to obtain an inclusive rate. For this reason, experimental cuts 
and the experimental infrared energy resolution can be implemented in our 
calculation, once the corresponding one-loop gluon Bremsstrahlung process is 
added to this calculation. 
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Abstract 

For Z — > 65, we calculate all the two-loop top dependent Feynman graphs, 
which have mixed QCD and electroweak contributions that are not factorizable. 
For evaluating the graphs, without resorting to a mass expansion, we apply 
a two-loop extension of the one-loop Passarino-Veltman reduction. This is 
an analytic-numerical method, which first converts all diagrams into a set of 
ten standard scalar functions, and then integrates them numerically over the 
remaining Feynman parameters, with rapid convergence and high accuracy. We 
discuss the treatment of infrared singularities within our methods. We do not 
resort to unitarity cuts of two-point functionsfor calculating decay rates; these 
are useful only to obtain an inclusive rate. For this reason, experimental cuts 
and the experimental infrared energy resolution can be implemented in our 
calculation, once the corresponding one-loop gluon Bremsstrahlung process is 
added to this calculation. 

High energy experimental data have reached such an accuracy that higher loop 
effects have to be accounted for. It is increasingly clear that new physics most likely 
lies beyond the rubric of tree or even one-loop graphs of the standard model. 

To asses these higher order effects, in quite a few processes where one can identify 
a single large scale, a very powerful technique is asymptotic momentum or large mass 
expansion 0, ^, This has been applied successfully in several instances, where 
either the external momenta are small compared to the internal masses, or vice versa. 
However, in many examples, the external momenta are comparable to the internal 
scales, and therefore an expansion in their ratios does not apply a priori, or at least 
is not economical for convergence. 

There is a further complication which must be addressed. Massless particles and 
particles with very small masses are with us. Partly to avoid potential infrared 
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and mass singularities, many authors calculate inclusive rates, where such issues are 
by-passed. On the other hand, in a typical experimental setup one needs to make 
kinematical cuts among other things, and therefore exclusive processes must be con- 
sidered. It would help if a theoretical calculation can accommodate this. 

In a previous article [0] we gave a general framework to calculate any two-loop 
amplitude with non-trivial interactions. Our approach fully respects the mass struc- 
ture and the kinematics of the physical process being investigated. The outcome is a 
standard set of ten functions and their derivatives, arising from tensorial decomposi- 
tion. Numerical integration over the remaining Feynman parameters is then used for 
obtaining final results. 

In this article we apply this general program to an 0{asg^) calculation of the 
process Z — ^ bb. So far, this process was calculated by mass expansion methods 
in refs. with the gluon Bremsstrahlung process integrated over the whole phase 
space. The calculations already available, being performed at higher order in the mass 
expansion, show that the expansion methods work well in this particular process. For 
this reason, it is to be expected that the numerical output of the calculation described 
in this paper will not result in a sizeable difference from the existing results of ref. 

Our emphasis here is more on an illustration, based on an important physical 
process, of how our two-loop methods work, because the whole machinery deployed 
covers almost all typical situations, including particles with various masses, two- and 
three-point functions, and treatment of infrared singularities. 

A complete two-loop calculation of the exclusive Z ^ bb decay consists of two 
technically different parts, namely the pure two- loop part, and the evaluation of the 
one-loop Bremsstrahlung process. In this paper we treat the mixed QCD-electroweak 
non-factorizable two-loop contributions to this process. Technically this is the most 
difficult part of a complete top-dependent 0{g^as) calculation because it is the part 
involving massive two-loop diagrams. Of course, for obtaining an infrared finite re- 
sult one needs to add the gluon emission graphs. The remaining factorizable and 
Bremsstrahlung diagrams can be calculated by conventional methods. The 0{g'^as) 
gluon emission process involves four-point one-loop diagrams, for which standard 
techniques such as the Passarino-Veltman reduction can be used. Here we concen- 
trate only on the pure two-loop part of the calculation, for which special massive 
two-loop methods are necessary. 

The coupling of the Z boson to a 6 quark pair is given by: 

^ n li^ivd-adi5) , (1) 
2 cos Ow 

which results into a. Z —>■ bb width given by the following expression: 

^z^bi= Mzi\v£ + \a,\') . (2) 

lOTT COS"^ tlw 

The tree level values a^j^^^'^^ = —1/2 and vj^''^^^ = — 1/2 + 2/3sin^6'vi/ receive both 
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QCD and electroweak radiative corrections. 

The two-loop graphs which contribute to the Z ^ bb decay to this order are shown 
in fig. 1. In the following we calculate these graphs for finite Z, W, and top masses, 
without using a mass or external momentum expansion. The only approximation 
we do in this calculation is to neglect the b quark mass, where this is justified. We 
use this approximation because its effect is very small, of the order of mh/Mz^w,t- 
This approximation simplifies the intermediary expressions resulting from the tensor 
reduction of the diagrams. However, we note that including a finite b mass throughout 
our calculation would be straightforward. 

The first step in calculating these diagrams, after some trivial Dirac algebra, is to 
decompose the resulting tensor integrals into a set of standard scalar integrals. This 
is done along the lines of ref . . All tensor structures can be reduced to a set of ten 
scalar functions 7ii-7iio, whose definitions are given in ref. 0. 

In carrying out this scalar decomposition procedure, there is a subtlety related to 
the apparent introduction of some spurious singularities. To illustrate this issue, let 
us consider a general tensor integral, which is of the following type (see ref. |^): 

{p'^ + m?)"i (g2 + [(r + + mi]"^ 

In the expression above, all momenta are rotated to Euclidian. Such a tensor integral 
can always be decomposed into Lorentz invariant integrals of the type: 

rq ivm<i-kf 

(p2 + mf )°i (g2 + m|)"2 [(r + kf + m|]"3 ^ ' 

by decomposing the loop integration momenta p and q into components parallel and 
orthogonal to the external momentum k of the two-loop integral: 

pl=p^-l^k^ , q^^ = q^^-1±k^ . (5) 

However, in doing this, one apparently introduces light-cone singularities of the 
type l/fc^. At the same time, it is obvious that the original tensor integral of eq. 3 
is free of such singularities, and in fact it can be evaluated analytically at vanishing 
external momentum for any combination of masses [0]. This obviously means that 
the scalar integrals which result from the tensor decomposition must have such a k"^ 
behaviour as to compensate the 1/k'^ singularities which result upon introducing the 
orthogonal loop momenta defined in eqns. 3. Indeed, one can convince oneself that 
1/P light -cone singularities are non-existent when the tensor integrals are reduced 
into scalar Tii functions. As an example, we give here the tensor decomposition of a 
two-loop integral with three indices: 

a pa q 



[{p + ky + mf]2 (g2 _)_ (j"2 _|_ jj^2\ 
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Figure 1: The three- and two-point two-loop Feynman graphs which contribute to the 
top- dependent correction of 0{asg'^) to Z bb. Diagrams Vi^s^g, 7,8,9 occur together 
with the diagrams obtained by reverting the fermion line, which are not shown explic- 
itly. 
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with the scalar functions Ti. and V defined in ref. [0]. By using the exphcit expressions 
of the 7i and V functions, one can easily see that this decomposition is indeed free 
of singularities at /c^ — > 0. Singularity-free decomposition formulae are obtained in a 
similar way for all other tensor integrals involved. 

We encoded all the necessary Dirac algebra into an algebraic manipulation pro- 
gram (we used both FORM and Schoonship). The computer program then uses the 
two- loop reduction algorithm which we described in detail in ref. 0. Thus it reduces 
the two-loop Feynman graphs into a set of standard scalar integrals. The final output 
of the algebraic program is a combination of the special functions hi, which are the 
finite parts of the two- loop scalar integrals Hi — see ref. 0. 

These final algebraic expressions are directly suitable for further numerical inte- 
gration. To do this, they are transfered into a FORTRAN integration package 0. 
The numerical package can calculate efficiently and with high numerical accuracy 
the functions hi, starting from their integral representations. Then it performs the 
numerical integration over the remaining Feynman parameters. For achieving both a 
high integration speed and a high accuracy of the final result, the complex singulari- 
ties of the integrand are found automatically, and then a smooth complex integration 
path is automatically calculated in terms of spline functions. Then, along this com- 
plex path, an adaptative deterministic algorithm is used, which leads to an accurate 
numerical evaluation of the Feynman graph. 

Collecting all proper vertex contributions shown in figure 1 (diagrams Vi-Vio), we 
obtain a correction to the Zbb vertex of the type 7^(1 — 75)^4'^^^''*^^) in the massless b 
limit. 

From the b self-energy diagrams we derive a b wave function renormalization con- 
tribution. Denoting the b self-energy by zS(p), we have in general: 



S(p) = A{p^) + 75^5(/) 1-P + B{p^) + -f^B^ip 



(7) 



For the general case of a finite b mass, the coefficients A(p^) and A^lp"^), and the mo- 
mentum derivatives of B{p'^) and B^{p'^) at p^ = ml are needed for extracting the b 
wave function renormalization contribution. In ref. we have shown how this can be 
done within our two-loop formalism. As a side remark, we note here that the on-shell 
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Figure 2: Extracting the infrared divergent pieces of the two-loop diagrams analytically. 
The infrared divergency of the two-loop diagram is the same as the infrared divergency 
of the product of the two one-loop diagrams obtained by "freezing" the common line 
in the loop momenta integration. 



external momentum differentiation can be best performed by working on the inte- 
grand before tlie momentum integrations 0. Tfiis reduces tlie size of the expressions 
which appear at intermediate stages of the calculation, and provides for a systematic 
treatment of all diagrams. The result for the two-loop wave function renormalization 
constant is again a set of the ten standard functions we have introduced. 

In the massless b quark limit, the expressions of the self-energy diagrams simplify 
considerably. The self-energy simplifies to = A(j9^)(l + 75)7 -p. Accordingly, the 
b wave function renormalization constant contribution to the Zbb vertex correction 

reads 7^(1 - 75){a'r + = 0). 

Therefore, the non-factorizable two-loop contributions of figure 1 contribute a 
correction 



7^(1 - 75) [A^""'''-'^ + {a'r + vT''^)A{p' = 0)] (8) 

to the Zbb vertex. 

At this point we would like to discuss the treatment of infrared divergencies. 
Among the three-point diagrams, V2, V4, Vg, and Vj are infrared divergent. Their 
infrared divergencies can be isolated analytically in the way shown in figure 2 for the 
case of V2. The analytical isolation of infrared singularities is based on the observation 
that the original two-loop diagram has one line common for the two loop momenta 
which are to be integrated out — the W/(f) line in the case of diagram V2. By 
"freezing" the loop integration momentum of this internal line to the value of the 
loop momentum whose one-loop sub-diagram is infrared finite — the ttW triangle 
sub-diagram in the case of diagram V2 — one obtains a product of two one-loop 
diagrams which has precisely the same infrared singularities as the two-loop diagram. 
All other infrared divergent diagrams can be treated in the same way. 
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The analytical infrared isolation and factorization into one-loop diagrams clearly 
works nicely for all two-loop diagrams involved in the Z ^ bb calculation, shown 
in figure 1. However, it is an open question if this can be done in all possible two- 
loop cases, and especially in the case of calculations of higher order in as, such as 
two-loop pure QCD calculations. This point clearly deserves further investigation. 
We note that our massive two-loop reduction scheme is mainly designed for massive 
calculations; for pure QCD calculations other methods which take advantage of the 
massless structure of the theory may prove more efficient — see for instance ref. |^ 
for a few recent examples. 

Once the infrared singularities are isolated and written in the form of one-loop 
diagrams, one chooses a regularization method to treat them and cancel them upon 
the real gluon emission process. Dimensional regularization is often used in higher- 
order QCD calculations. Given that the process considered here is only of 0{as) in 
the strong coupling constant, a second possibility is to use a gluon mass regulator. 
This does not upset the Slavnov- Taylor identities for our case because to the order 
considered here the infrared structure is the same as in the Abelian case. 

Independently of the choice for the infrared regulator, separating the infrared 
singular piece analytically is useful for increasing the efficiency of the numerical in- 
tegration of the two-loop diagram. For instance, in the case of diagram V2 discussed 
in figure 2, the infrared singularity appears as an end-point singularity in a two- fold 
Feynman parameter integration. By subtracting first the product of two one-loop 
diagrams shown in figure 2, this end-point singularity is being removed and the nu- 
merical integration over the resulting smooth function (V2 — ^2*^^^"*) becomes much 
more efficient than V2 alone. 

In the following discussion it is understood that the infrared singularities are 
regularized by introducing a gluon mass regulator nig. Again, we would like to stress 
that dimensional regularization can be used equally well. 

When calculating the two-loop diagrams of figure 1, we deal with infrared diver- 
gences in the form of Lirrimg^o ln{mg), where rrig is the gluon mass regulator. These 
infrared divergent terms cancel those from real gluon emission when we calculate a 
decay rate, and we have checked this explicitly. The net effect of this cancellation is 
to replace the gluon mass regulator rrig by the maximum undetected gluon energy, or 
the energy resolution of an experiment uJmax — the logarithmic infrared behaviour is 
universal and we do not need to calculate explicitly the Bremsstrahlung process to 
extract it. Of course, after the infrared cancellation one is left with a finite piece and 
an explicit calculation of the real gluon emission process is needed. This involves the 
evaluation of four-point one-loop diagrams and integration over three-particle phase 
space with the appropriate experimental cuts. The Bremsstrahlung process is beyond 
the scope of this article — we only note that the calculation can be performed with 
conventional methods such as the Passarino-Veltman reduction. 

In figure 3 we give numerical results for the pure two-loop contribution to the 
Z ^ bb amplitude, stemming from the diagrams shown in figure 1. Instead of ploting 
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Figure 3: Numerical results for the pure two-loop contribution due to the proper Feyn- 
man graphs given in figure 1 (see the text). The ultraviolet poles are removed by min- 
imal subtraction. The infrared divergencies are separated as one-loop integrals and 
can be handled either by dimensional regularization or by using a gluon mass regu- 
lator rrig. Here we give numerical results for a range of nig. An overall colour and 
coupling constant factor of ias{g^ /V2 cos 9 w) is understood. 



each diagram separately, we give their total contribution, which translate into the 
term _|_ {^ol'p^ + v'^^'^^^')A{p^ = 0) of eq. 8. All diagrams are subtracted in 

the ultraviolet by minimal subtraction. We have checked explicitly that the ultraviolet 
infinities are absorbed into the mass, coupling constant, and wave function renormal- 
ization. This serves as a good check on the two-loop tensor decomposition algebra. 
The infrared divergencies are regulated by using a gluon mass regulator. Should di- 
mensional regularization of infrared singularities be preferred, this is straightforward 
to implement by correspondingly treating the one-loop infrared pieces of the type 
shown in figure 2. We used Mz = 91.187 GeV, Mw = 80.41 GeV, and the effective 
electroweak mixing angle sin^ 9 = .2312. 

In conclusion, we calculated all two-loop proper diagrams relevant for the 0{asg'^) 
top-dependent correction to the Z ^ bb decay width, by using two-loop methods 
which we developed previously. We have developed an on-shell momentum expansion 
which allows for a simple extraction of wave function renormalization constants. We 
discussed the isolation of infrared singularities in the form of products of one-loop 
integrals. This allows the treatment of infrared cancellations with the Bremsstrahlung 
process either by dimensional regularization or by using a gluon mass regulator which 
is legitimate in this particular order. Our calculation is exact, in the sense that we 
do not rely upon mass or momentum expansions of Feynman diagrams — we rather 
calculate them at finite mass and momentum values from the outset. Being based 
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on a general algorithm, the calculation presented in this contribution opens the way 
for the calculation of a number of other two-loop corrections which have become 
necessary because of improved experimental electroweak data. 
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